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Abstract 

We investigate, within the framework of three generations of neutrinos, the 
effects of CP violation in long-baseline neutrino oscillation experiments. We 
aim at illuminating the global feature of the interplay between genuine effect 
due to the CP violating phase and a fake one due to the earth matter effect. 
To this goal, we develop a formalism based on the adiabatic approximation 
and perturbative treatment of the matter effect which allows us to obtain 
approximate analytic expressions of the oscillation probabilities. We present 



an order-of-magnitude estimation and a detailed numerical computation of 
the absolute and the relative magnitudes of the CP violations under the mass 
hierarchy suggested by the atmospheric neutrino anomaly and the cosmolog- 
ical dark matter. We find that the genuine CP violating effect is at most ~ 
1 %, and the matter effect dominates over the intrinsic CP violation only in 
a region of parameters where the oscillation probability of —>■ v^. is large. 
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I. INTRODUCTION 



The origin of CP violation is still a mystery in particle physics. Unlike in the quark sector 
even the very existence oi CP violation is not known in the lepton sector. Recent advances 
in neutrino observations mainly of astrophysical origins strongly suggest the existence of 
tiny neutrino masses [|T],0. If this is the case nature would admit the Cabibbo-Kobayashi- 
Maskawa (CKM) |^ type flavor mixing also in the lepton sector which allows us to have 
CP violating phases with three (or more) generation of leptons. If revealed, it should give 
us important insight into our understanding of fundamental structure of matter. Moreover, 
there is an intriguing suggestion @] that CP violation in the lepton sector might be an 
indispensable ingredient in producing the baryon asymmetry in the universe. 

It has been known since long time ago that the effect of CP violating phase can 
in principle be observable in neutrino oscillation experiments. The other earlier references 
include As we will recapitulate below the particle-antiparticle difference between 

the oscillation probabilities, /S.P{uj3 Ua) = P{i^f3 — ^ i^a) — P{^i3 ^a), in vacuum is 
characterized by the leptonic analogue of the Jarlskog factor, the unique and phase- 
convention independent measure for CP violation 0. 

Recently, measuring leptonic CP violation in the long-baseline neutrino oscillation ex- 
periments |TU|-[T^ received considerable amount of attention in the literature The 
potential obstacle in measuring CP violation in the long-baseline experiment is the matter 
effect in the earth. It is well known that the matter effect acts differently in propagation 
of u and u in matter; it gives rise to the index of refraction which differs in sign between 
u and u HTSf. Then, the u — u difference AP of the oscillation probabilities is inevitably 



contaminated by the matter effect |]I6[. In fact, it is known by numerical computation that 
the matter effect is overwhelming over the genuine CP violating effect at certain values of 
the mixing parameters in the z/^ Ue channel But it appears to the authors that we 

still lack understanding of the over-all features of the relationship between the CP violations 
due to the matter and to the CP violating phases. 
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It is the purpose of the present paper to illuminate global structure of the interplay 
between the matter and the genuine CP violating effects in long-baseline experiments. To 
this goal we develop a formalism by which we can derive approximate analytic expressions 
of oscillation probabilities. These analytic formulas will allow us to have global view of the 
features oi CP violation in neutrino oscillation experiments. Our formalism is based on the 
adiabatic approximation and takes into account the matter effect in a perturbative way. It 
also enjoys further simplification due to the presumed hierarchy of neutrino masses which 
will be explained below. 

This is the first in a series of papers to investigate the same problem in various neutrino 
mass hierarchies. In this paper, we assume that the atmospheric neutrino anomaly, which 
is observed by the Kamiokande, 1MB, and the Soudan 2 experiments 0], can be interpreted 
as the evidence for neutrino oscillations whose relevant mass scale is Am? ~ 10^^ eV^ or 
larger. (We note, however, that the two experiments do not observe the atmospheric neutrino 
anomaly [1^.) It is a natural assumption because it is the very motivation for planning 
the long-baseline experiments. The restriction leads to the hierarchy between the matter 
potential and the mass differences, which allows us to treat the matter effect perturbatively. 
Namely, we derive closed-form analytic expressions of the neutrino-antineutrino difference 
between oscillation probabilities which is generally valid under the adiabatic approximation 
and a first-order perturbative treatment of the matter effect. 

We restrict ourselves into the case of neutrino mass hierarchy motivated by the cosmo- 



logical hot dark matter |jT8[ in our analysis of the features of CP violation in this paper. 
The assumption of the dark matter scale mass difference allows us to utilize the strong con- 
straints on the mixing parameters deduced from the terrestrial experiments It is the 
reason why we decided to investigate this case first; Because of the constraints we are able 
to estimate the magnitude of CP violation in a less ambiguous way. 

We give a few ward on the problem if the mass hierarchy we examine in this paper can 
account for the other known hints for the neutrino masses. They include the solar neutrino 
deficit and the LSND experiments |20|- Let us first discuss the solar neutrinos. It is well 



known that for this scheme the capabihty of explaining the gross deficit of the solar neutrino 
flux depends upon the parameter regions. In the region (A) (see Eq. (^) for the definitions 
of these regions) the global deficit of the flux can be accommodated in a manner of Acker 
and Pakvasa I^Tj. In the region (B) there is no way of having a gross deficit of the solar 
neutrino flux within the framework of three-flavor mixing. Therefore, we must introduce 
some other ingredients such as sterile neutrinos to achieve the deficit in the region (B). 

The results of the LSND experiments can be easily accommodated in the present frame- 
work. The LSND data allow the mass scales which is consistent with hot dark matter. If 
smaller values of Am^ is preferred for consistency with other experiments it may be necessary 



to relax the constraints on the mass scale required by the mixed dark matter cosmology [jTH 
If the simultaneous consistency with the solar neutrino and the LSND data is demanded 
there is no solution within our present framework. 

Within the restrictions to the mass hierarchies hinted by the dark matter and the atmo- 
spheric neutrinos, we will try to answer the following questions: 

(1) Are there any channels which are much less contaminated by the matter effect? 

(2) In which parameter regions do we expect to have maximal CP violation and how large 
is its size? 

It is known that the matter effect in the oscillation probability in long-baseline exper- 
iments is not very large, at most a few to several %. One might then feel strange that 
the matter effect is dominant in certain channel. The point is that what we are dealing 
with is not the oscillation probability itself but the difference between the neutrino and the 
antineutrino oscillation probabilities. The matter effect can give rise to a dominant effect in 
such V — v difference. 

This investigation was motivated by a question asked by a long-baseline neutrino ex- 
perimentalist [Q; "Is the matter effect contamination small in channel?" This is 
the interesting question for two reasons, one theoretical and one experimental. Experimen- 
tally it is a very relevant question because they are planning to do long-baseline neutrino 
experiments in the appearance channel, i/^ — z/^. It is also of interest from the theoretical 



point of view; In the conventional treatment of "optional" two-flavor mixing favored by ex- 
perimentalists the Vfj, Vt channel might be free from the matter effect because there is 
no /x's and r's in the earth. On the other hand, CP violation is the genuine three-flavor 
mixing effect which cannot occur in two-flavor mixing framework. Therefore, the question 
still remains whether the matter effect is small in the — > 1/^- channel, as correctly raised 
by the experimentalist. 

In Sec. II, we review some basic facts on CP violation in vacuum in the context of 
neutrino oscillation experiments. We also summarize the mass patterns of neutrinos which 
we use in our analysis in this paper. In Sec. Ill, we set up our formalism based on the 
adiabatic approximation. In Sec. IV, we develop a framework for perturbative treatment 
of the matter effect which is applicable to the long-baseline neutrino experiments. In Sec. 
V, we derive the approximate analytic formulas for the neutrino-antineutrino difference in 
oscillation probabilities by taking account of the matter effect to first-order in perturbation 
theory. In Sec. VI, we use our analytic formula to illuminate global features of the competing 
two effects producing CP violation, the genuine effect due to CP violating phase and a fake 
one due to the matter effect. In Sec. VII, we give the results of our detailed numerical 
computation of the CP violation to confirm the qualitative understanding of the structure 
of the coexisting two effects obtained in Sec. VI. The last section VIII is devoted to the 
conclusion. In Appendix we confirm that the adiabatic approximation is in fact a very good 
approximation for matter density profiles relevant for the long-baseline experiments. 

II. CP VIOLATION IN VACUUM AND NEUTRINO MASS SPECTRUM 
MOTIVATED BY DARK MATTER AND ATMOSPHERIC NEUTRINO 

ANOMALY 

In this section we first review briefiy the CP violation in neutrino oscillations in vacuum. 
We then describe the neutrino mass spectrum which we consider in this paper. It is intended 
for convenience for the readers who are not familiar to the subject and precedes the following 
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three sections in which we develop our analytic framework usable for analyses with more 
general neutrino mass patterns. 

We work with the three-flavor mixing scheme of neutrinos and introduce the flavor mixing 
matrix U as 



(1) 



where z/Q,(a = e,fi,T) and z/j(z = 1,2,3) stand for the gauge and the mass eigenstates, 
respectively. Then in vacuum the direct measure of CP violation can be written 0,^] as 



AP{u0 Pa) = P{^p T^a) - P{yfi = -4J/3Q(sin Ai2 + sin A23 + sin A 



31 y 



(2) 



where 



A,; 



{rrij - mf)L 
2E ' 



(3) 



with mi{i = 1,2, 3) being the mass of i-th neutrino, E the neutrino energy and L the distance 
from the neutrino source. We have used in @ the Jarlskog parameter 



J/3a = Im[[4i[/*2f/^if//32] 

which is unique, up to the sign, in three-flavor neutrinos. 

We will use the following form of the neutrino mixing matrix 



(4) 



U 



C12C13 



'S12C13 



AS 



"'S12C23 — Ci2S23Si3e*'' C12C23 — Si2S23Si3e*'^ S23C13 



45 



AS 



(5) 



S12S23 - Ci2C23Si3e'" -C12S23 - 5120235136"" C23C13 

which is identical with the standard CKM matrix for quarks. Here Cjj = cos 9ij and Stj = 
sinOij where = (1,2), (2,3) and (1,3). With this parametrization Jg„ can be expressed 
as 



Jf3a = ±J, J = Ci2Si2C23S23Ci3Si3Sin(5, 



(6) 



where -|- sign is for cychc permutations, i.e., (a,/?) = (e, yu), (yU, r), (r, e) and — is for anti- 
cyclic ones. 



Let us come to the mass hierarchy which will be used in our analysis presented in Sees. 
VI and VII. We work with the mass hierarchy 

lAm^l ~ IQ-W^ Am^ = - ml, (7) 
lAM^I > a few eV^ AM^ = - - - m^, (8) 

where the first is suggested by the atmospheric neutrino anomaly [§, and the second is 
motivated by the hot and the cold dark matter cosmology 

The Kamiokande, 1MB, and the Soudan 2 experiments observed a 30-40 % deficit in the 
double ratio {^)observed/ expected [0- A natural interpretation of the anomaly is due to the 
neutrino oscillation. In particular, the multi-GeV data from the Kamiokande experiments 
indicates that the double ratio has a zenith-angle dependence which is quite consistent with 
the interpretation of the anomaly as the evidence for neutrino oscillations with the mass 
difference Am^ ~ 10"^ eV^. 

The hot and the cold dark matter cosmology is one of the viable models of the structure 
formation in the universe |TB[. The hot dark matter is the indispensable ingredient in the 
scenario by which the magnitude of density fluctuation normalized by COBE data can be 
made consistent with that of small scales determined by correlations between galaxies and 
clusters . Neutrinos of masses 2-20 eV are the natural candidate for the hot dark matter. 



In fact, it is the only known particles among the candidates for particle dark matter. 

The advantage of the dark matter motivated mass hierarchy (H) is that the mixing 
parameters are subject to the powerful constraints that comes from the reactor and the 
accelerator experiments [0] and we can draw a clear answer to the question on how large 



is the magnitude oi CP violation. 

Under the assumption (|) the oscillations due to the larger mass squared difference AM^ 
are averaged out in long-baseline experiments and we have, 

/ Am^L\ 

AP(z/^ ^ z/„) ~ -4 J,^ sin I J (9) 

Since (^) = 2.5 (^^^) (^) (i^)"' can be of order unity in long-baseline 



experiments, the CP violation APf^a can be approximated as ~ —AJafs- 

Let us estimate how large CP violation can be for our choice of mass hierarchy (|^) and 
In Fig. 1 we present equal- J contours on tan^ ^^13 — tan^ ^23 plane with use of the 
parameters 612 = it/ 4: and 6 = n/2 which maximize J. In the same figure we also plot 
the regions of parameters excluded by the reactor and accelerator experiments for the case 
|AM^| = 5 eV^ obtained in Ref. [^]. We notice from Fig. 1 that allowed parameters are 



restricted into three separate regions. 

(A) small — Si3 and small — S23 

(B) large — S13 and arbitrary — S23 (10) 

(C) small — Si3 and large — S23 

The constraints from the reactor and accelerator experiments are so strong that one neutrino 
flavor almost decouples with the remaining two. Namely, the Ur, Ve and almost decouple 
with the other two flavors in the regions (A), (B) and (C), respectively. The remaining two 
neutrinos can be strongly mixed with each other in each allowed region. We observe that 
the value of 4 J is at most ~ 0.01 in the regions (A) and (B) and ~ 0.001 in the region (C). 

In Fig. 2 we show the neutrino mass spectrum which is realized in the region (A) and 
(B). We do not consider in this paper the region (C) since the atmospheric neutrino anomaly 
cannot be accounted for and moreover CP violation is very small in this region. There exist 
two different mass patterns which can be realized in each parameter regions. Depending 
upon the sign of AM^ the decoupled state z/3 can be the heaviest, (A-1) or (B-1), or the 
lightest, (A-2) or (B-2). The case (A-1) is theoretically appealing because of the seesaw 
mechanism pH] and from the observed mass hierarchy among charged leptons. 



III. ADIABATIC NEUTRINO EVOLUTION IN MATTER 



We discuss the neutrino propagation in the earth matter within the framework of three- 
flavor mixing of neutrinos. It is worth to note that when we discuss the long-baseline 
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experiments whose baseline distances less than 1000 km the adiabatic approximation is a 
very good approximation, as we demonstrate in Appendix. It is because the neutrino beam 
only passes through thin layer of the continental structure closest to the earth surface in 
which the matter density is approximately constant, ~ 2.72 g/cm^ ||2^. By virtue of this 
fact we can derive a general closed form expression of AP, the i> — u difference in oscillation 
probabilities. 

The evolution equation of neutrinos (and antineutrinos) can be written in terms of the 
gauge eigenstate as 



. d 

dx 



where the Hamiltonian H{x) is given by 











= H{x) 











(11) 





ml/2E 










a(x) 








H = U 





ml/2E 





f/+ + 



















ml/2E 














(12) 



where x is the position along the neutrino trajectory. Here a{x) represents the matter effect 
and has the form 



a[x] 



±V2GfN, 



(13) 



where Gp is the Fermi constant, Ne{x) is the electron number density at x in the earth, and 
the + and — signs are for neutrinos and antineutrinos, respectively. 

We introduce the unitary transformation V{x) which diagonalize H[x) locally, i.e., at 
each point x of neutrino trajectory; 



V+{x)H{x)V{x) = Hd{x) 



(14) 



and parametrize the diagonalized Hamiltonian as Hd{x) = diag.[hi{x), h2{x), h^i^x)]. We 
then define the matter-mass eigenstate, mass eigenstate basis in matter, as 
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= V+{x) 











(15) 



The neutrino evolution equation in terms of the matter-mass eigenstate takes the form 



.d_ 

dx 



+ i^V+{x) -Vix) 
dx 



■'7712 



■'7713 



(16) 



The adiabatic approximation amounts to ignore the second term in the right-hand-side of 
eq. (p!6|). We will return in Appendix to the question if it gives a really good approximation. 
Under the adiabatic approximation it is straightforward to solve the evolution equation 

as 

(17) 



Va^x) = Vai{x)exp <{ -i dx'hi{x') ^ V^* (0)z/^(0). 



We then obtain the probability of the neutrino oscillation z/^ Va where vp is created at 
a; = and Va is detected aX x = L. It reads 



id 

xexp J dx[hi{x) — hj{x)] 



(18) 



We assume in this paper an idealized situation where the matter densities at the pro- 
duction and the detection points are identical. This should give a good approximation to 
the real experimental situation because these points are either on or are very close to the 
earth surface. One way argue that one can take V^(0) = V{L) equal to the vacuum mixing 
matrix U by saying that we design the experiment so that the production and the detection 
points of neutrinos on the earth surface, namely, in vacuum. We argue that even with such 
experimental condition it is the better approximation to take as V{0) = V{L) the value 
of matter-mass mixing matrix with matter density ~ 2.7 g/cm^ at just below the earth 
surface. First of all, if we take the vacuum mixing matrix U for V{0) = V{L) we have to 
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worry about the failure of the adiabaticity condition at the earth surface. The oscillation 
length of neutrinos is approximately given by 

— (iciv)!^)"'- («) 

and is much longer than the decay tunnel or the detector hole for interesting regions of 
Am^. Therefore, the neutrinos do not know if they take off at the earth surface prior to 
the detection and would rather feel as if they remain in the earth matter at the point of 
detection. 

By separating the summation into i = j and i ^ j terms one can rewrite the expression 
of oscillation probability into the form 

P{uf, ^u^;6,a)^ -45:i?e[K,(a)TC^.(a)y;,(a)y^,(a)]sin2[^7,,(a)] 

j>i 

+2J2MVa^{a)V:^{a)V^\{a)V^,{a)]sm%,{a)], (20) 



where 



Iij{a) = / dx[hi{x) — hj{x)]. (21) 

t/ 



IV. PERTURBATIVE TREATMENT OF MATTER EFFECT 

Now we intend to evaluate AP to first-order in matter perturbation theory. Let us 
first confirm that the perturbative treatment is reliable under the mass scale with which 
we are working. Since we are assuming the mass difference Am^ ~ 10~^ eV^ relevant for 
atmospheric neutrino anomaly we have, 

^ - 10-^ ( ^ (^] " eV (22) 

On the other hand, we estimate the matter potential for the continental structure as, 

a = 1.04 X 10-'' ( ^] ( ^) eV, (23) 

V2.7g/cmV VO.Sy ' ^ ^ 
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where Yg = Np/ {Np + A^„) is the electron fraction. Hence we see that the hierarchy between 
the energy scales, a -C (atmospheric mass scale), holds. 

We use stationary state perturbation theory in accord with the slow variation of a{x) 
on which the adiabatic approximation is based. It is convenient to take the vacuum mass 
eigenstate Uili = 1,2,3), which diagonalizes the unperturbed Hamiltonian, as the basis of 
matter perturbation theory. We denote the first and the second terms of the Hamiltonian 
(|T^) as Ho and H', respectively. We use the vacuum mass eigenstate Ui defined by 

i/,(x) = e-"''°'^u„ (24) 

where hf'^ = mf/2E, as the basis of matter perturbation theory. 

The Hamiltonians in this basis are Hq = U^HqU and H' = U^H'U, respectively. Since 

{H')a/3 = a6aeSf3e WC haVC 

{H')j, = aU,,U:^ (25) 
Then, the matter mass eigenstate can be expressed to first order in a as 

'mi - '^i + Yl , (0) A<d)^i (2^) 



It can be converted to the expression of vacuum mass eigenstate expressed by the matter 
mass eigenstate 

E^UgjUej fo'7\ 
TTq) 7(0)'^™-' ^ ^ 

j^i lli — 11 j 



Since the flavor eigenstate has the simple relationship with the vacuum mass eigenstate as 
= Uaii'i we obtain 



E 



TJ ■TJ*.TJ ■ 
^aiOij (0) _ (0)*^ 



(28) 



The matrix element of V can be read off from this equation as 

13 



Kfi — Uai ^ 7,(0)'^ 



It is instructive to verify that Vai in (^) satisfies the unitarity to first order in a, as it 
should. 

We have obtained the matrix V in the form V = U + 6V, where 6V denotes the correction 
first order in a. Then we can write, schematically, WW in (^) as 



WW = UUUU + UUU5V 



(30) 



where the second term actually contains four pieces. It is important to know the symmetry 
property of these terms. It follows that 



Re{UUUU) 
lm{UUUU) 

MUUU6V) ^-zT^r^ 



lm{UUU5V) 



Re{UUU5V) 
-lm{UUUU) 
-Re{UUU6V) 
lm{UUU6V). 



They stem from the fact that UUU6V is linear in a and that the transformation 6 
equivalent to U ^ U* . 



(31) 



-6 is 



V. CP VIOLATION EFFECT IN THE PRESENCE OF MATTER 



The neutrino-antineutrino difference can be written as 



= P{ui3 ^ z/q ; 5, a) - P{ui3 fa ; -5, -a) 



(32) 



Thanks to the symmetry properties (pT]), we can express eq. (|32D as 



AP(z/^ ^ z/„) = - 4^Re(f/f/f/f/) 



+ 2Y,\m{UUUU) 



al3 ; ij 



j>i 



sin^{-/ij(a)} - sin^ {-Iij{-a)} 



sin lij (a) + sin Iij{—a) 



-4:J2^e{UUU5V)^i3. 



sin {-Iij{a)} + sin {-hji-a)} 



— sin/ij(— a) 



(33) 
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where 



UaiU*jUpiUi3j 



(34) 



etc. We note that the last term in (|3^) is at least second order in a. 
The energy eigenvalue hi can also be obtained to first order in a as 

2 

hi = -^ + \Uei\'^a. 

Therefore, Iij{a) defined in (|2l| ) can be given by 

hM) = + (|f/«p - |f/e,f ) / dxa{ 

zHj Jo 



X) 



(35) 



(36) 



where Am?- 



mf. Up to this point one discussion relies only on the hierarchy 



a -C \Amfj\/E and the hierarchy among Amfj need not to be assumed. 

To obtain an explicit form of AP we need the expression of UUU6V. Due to the mass 
hierarchy and (H), we can combine the terms and calculate J2i=i,2{UUU6V)a,i3 -is by 
ignoring the difference between Am^g and Am^^. We also calculate the sum {UUU5V)a^i3- 12 
separately: 

AEa 



J2iUUU6VU.,,3 

i=l,2 



AM2 



u, 



e3| 



2|?7a3nf^/33p 



a3\ 



(37) 



{UUU5V)^p,,2 
2Ea 



Avrfi 
2Ea 



U:,Ue2UalU:,{\UfS2\' - \U^in + UelU:^U*M2{\Ua2r - \Ual 



AM2 



{31 



In these equations we have dropped the terms further down by -^f^- Notice that the terms 
containing Am^ cancel out in J2i=i,2{UUU6V)ai3 -is- Note also that only the mixing matrix 
elements Ua3 appear in ( pTf ) as it occurs in the oscillation probability in vacuum. 

Combining all these together we obtain the expression of AP which is valid to first order 



m a: 
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— 4Jcos 



f/elP — \Ue2\^ ] / dxa{x) 



smi L 



+ARe{U^^Ul^U;^Up2) sin 
-Re 



+t/|lf//32f/elt/:2(l^:^a2|' 



f^alf^:2f^:if^e2(|f//32|'-|?7^l| 



2 



Am^ ^ 



sm 



Am^ 



AM2 



Re 



32Ea 



/33'^el^e3 
2 



e3^ 



sm 



Am^ 



-L 



f/esr 



sm 



,/AM2 



(39) 



where the term proportional to lm{UUU6V) and the a-dependent piece in the last term in 
(39) are ignored because it is of order or higher. 

A few remarks are in order concerning the sign of AM^ and Am^. Strictly speaking, 
the matter effect distinguishes between the neutrino mass spectra of the types (A-1) and 
(A-2), or (B-1) and (B-2) defined in Fig. 2. However, with our choice of mass hierarchy, 
AP^up — > Ua) barely depends on AM^, and consequently the final results do not depend on 
the sign of AM^. It is because the last two terms in (^) can be safely neglected due to the 
hierarchy in energy scales, 

-1 



Ea 



1.04 X 10" 



P 



E \ /AM2^ 
IGevJ lleV^ 



(40) 



AM2 ■ V2.72gcm-3^ 

Hence, we do not distinguish between the mass patterns (A-1) and (A-2), or (B-1) and (B-2) 
in this work. (Note, however, that they can be distinguished by consideration of r-process 



nucleosynthesis in supernova 



For the smaller mass difference, Am^ 



TTln 



mf, we 



assume that it is positive. But, it is easy to accommodate the case with negative Am^ in 
our analysis after ignoring the last two terms in (^), which we will do in our subsequent 
analysis. All we have to do is to change the over-all sign of AP(z/^ — > //„). 

By using the explicit form of the mixing matrix (^ it is then straightforward to obtain 
AP, the neutrino-antineutrino difference between oscillation probabilities, for — > and 



i^fi channels: 
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P(iy„ 



P(l7„ 



— 4Jcos 



rL 

c^3 COS 26^12 / dxa{x) 
Jo 



sm 



2E 



+ 
+4 



- sin^ 26*1 ^c?, f Co 
Ea 



_ , _ _ 2 2 
12Ci3l,C23 ■S23'^13 



C^3 COS 2^^ 



4 cos 2^12 J cot 5 sm 
sin 26^12 sm46'i2C?3(c23 — S23S?3) + 4cos46'i2Jcot5 



19 / dxa{x) 



sm' 



sm 



Am^L 
2E 



AE 



(41) 



AP(z/, 



+ 



4 J COS 

1 + ^2 



Ci3 COS 26* 



12 / dxa{x) 





/Am'^L 



2E 



sin|^- 

^ cos 26*12 cos 26*23 - 4:Cy^J cot 6}^ 



+ ;cf3(sin2 2^12 + sin^ 2^23) - 



l + .?3^^• 



sm 



Ci3 cos 2^12 / dxa{x] 
Jo 



^ Ea 
Am^ 



sin 26*12013 {sin 4^i2 sin 26*23 



2 ofl / 1 + ■^is 



2 \ 2 



- sin4^i2Si3} 



-4 cos 46*12 cos 26* 



^23 



1 + s 



2 

V cot (5 + 32 cos 2^12^3^^ cot^ S 



sm 



Am^L 
2E 



sm 



Am^L 
4E 



For completeness we also give the expression of AP for Ue Vt channel: 



(42) 



2 /I /■'^ 7 / n1 f Ani?L 

Ci3 COS 26*12 / dxa[x 
Jo 



— 4Jcos 



+ 

+4 



/o 

- sin^ 26*i2Ci3(s23 — C23S13) + 4 cos 26*12 J cot 
Ea 



sm 



V 2E 

12^ 



sm 



Ci3 cos 26* 



12 



sin 26*12 sin4^i2C?,(soq — Co,s?,) — 4 cos 46*1 2 J cot 5 



/12<-i3l'^23 



23 ■'13^ 



'12^ 



sm 



(ixa(x) 

Arri^L 
AE 



sm 



Am^L 
2E 



(43) 



Let us compare these analytic results with the exact solutions obtained in pSf for a 
constant matter density to examine the accuracy of our approximate formulas. We pick up 
the following two parameter sets, 



a s 
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3.0 X 10"^s^3 



(b) 5^3 = 2.0 X 10-2,4 



2.0 X 10^2 
0.98 



(44) 
(45) 



from the allowed regions (A) and (B), respectively. These parameters are chosen so that J 
takes a maximal value within the each allowed region and are plotted in Fig. 1. We will use 
the same parameter sets also in our analyses in the following sections. To obtain the "exact 
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results" by using the analytic expressions given in ||28| we take an average over the rapid 
oscillations due to the larger mass difference AM^. We take the constant matter density 
and electron fraction, p = 2.72 g/cm^ and Yg = 0.5. In Fig. 3 we show the comparison 
between the exact and approximated values of AP(z/^ —>■ z/g) and AP(z/^ —>■ Ur). We have 
fixed the remaining parameters as — 0-3; AM^ = 5 eV^ and 6 = 7t/2. We see that for 
~ 1 GeV the approximation is very good even for smaller values of Am^ (~ lO^^eV^ ) 
than the one we assumed in eq. (0). 

VI. CP VIOLATION VS. MATTER EFFECT IN LONG-BASELINE NEUTRINO 
OSCILLATION EXPERIMENTS: AN ORDER-OF-MAGNITUDE ANALYSIS 

In this section we present the results of our order estimations of the magnitude of CP 
violation, and give an answer to the question of relative importance between the genuine 
CP-violating and the matter effects. This precedes the presentation of the results of detailed 
numerical computation in the next section, which will provide us complementary informa- 
tions. The order-of-magnitude estimation based on our analytic expressions of AP illumi- 
nates the general features of the interplay between the CP phase and the matter effects, and 
is valid in whole allowed parameter regions. On the other hand, the numerical computation 
using the exact expression of the oscillation probability (albeit under the constant density 
ansatz) will reveal precise features of the relative importance between competing two effects. 

For convenience of our discussion let us denote the three terms in (^l]) and (^) as 

AP{u(s ^ z/J = APgP + APpf'^"^' + APpf''^'^ 

^ Cgf sin A + C^f'^^' sin A + C^f'^^'' sin' (A/2) , {f3 = p, a = e,T ), (46) 

where A = Am^L/2i?, which can be of order 1 for Am^ < 10^^ eV^, E >1 GeV and 
L > 100 km. C^ is the coefficient of the first term in these equations which represents 
the genuine CP violating effect corrected by the matter effect. C^'***^''^ and C^'***^'^^ are 
the coefficients in the second and third terms in (^TJ) and (0), display the matter effects 
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from the different sources. That is, C^'^**^''^ represents the matter effect which arises due 
to the evolution of the phase of the neutrino wave function in matter, whereas C*^"**'^''^ 
comes from the correction to the mixing matrix, as exhibited in the UUU6V factors in (p3[). 
We notice that C^^ and CJ^^^^'^^^ contain the integral Jq dxa{x) which can be estimated, 
by assuming the constant density 2.72 g/cm^, as aL = 1.2 x 10"^ for KEK-PS— >Super- 
Kamiokande experiment with baseline of 250 km, and aL = 3.5 x 10^^ for MINOS or 
CERN-ICARUS experiment with baseline of 730 km or 732 km, respectively. On the other 
hand, C^""^''^ carries the coefficient Ea/Am'^ which is of the order of ~ 10"^ for E ~ 1 GeV 
and Am^ ~ 10~^ eV^. Therefore, apart from prefactors which depend on mixing angles, we 
estimate that 

Cgf oc -4J ~ 10-^ (47) 
C™'^""'^^ oc aL ~ 0.1, (48) 



We recognize that C^"**^'"^ can be much larger, depending on the mixing angles, than the 
genuine CP violating effect C^^. We also note that from eq. ( ^Qf ) it is clear that if we 
assume smaller Am? values C^"**^'"^ become larger (See also the Tables 1 and 2). 

In Table 1 we summarize the results of our order-of-magnitude estimation of these three 
coefficients for z/^ — > v^. channel. In doing the estimation we have taken into account the 
coefficient which depend on the mixing angles, and the numbers presented in Table 1 refers 
to the possible maximal values in each region. We notice that, if sin A ~ 1, in the region 
(A) the oscillation probability -P(t'/i — > Ve) can be large, ~ 1, but CP violation due to the 
matter effect is much larger than the intrinsic CP violation effect, AP^^ ~ 0.1 x /\p»T*ai*er_ 
On the other hand, in the region (B), P{Vfj, Ve) is small, ~ 10~^, but the contamination 
of matter effect in CP violation is very small, AP'"'***^'^ ~ 10~^ compared to AP^^ ~ 10~^. 

In Table 2, we present the same quantities for channel. We observe that in this 

channel the intrinsic CP violating effect is larger than the matter effect in both regions (A) 
and (B), i.e., AP^^ ~ 10"^ whereas AP™<^**^^ ~ lO'^. 
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Let us try to understand the qualitative features of these results. Probably, the most 
interesting aspect of our results is that the matter effect is small in the region (B) in z/^ — > Vf. 
channel. But, in fact it is not difficult to understand the reason why. In the region (B) 
the mixing parameters are such that v^. is much heavier than z/^ and z/^, which are almost 
degenerate, or that and are heavier than v^. by the same amount in the mass hierarchies 
given in Fig. 2. The matter effect only affects electron neutrinos and the matter potential 
is small compared with Vg-Vy, mass difference, IS.M'^ jE ^ a. Then, it is easy to expect that 
the matter effect is small not only in —* but also in Vf, channels. 

In the region (A), on the other hand, is much heavier than Vf, and z/^, or vice versa, 
and z/e and are strongly mixed. Therefore, one would naively expect that the matter 
effect is sizable, as it is indeed the case in z/^ v^, channel. But, the situation is different 
in Vt- channel. A z/^ can easily communicate with z/g and thus feels the effect of earth 

matter but to make oscillation into z/^ it has to overcome the huge mass difference compared 
with the matter potential. Therefore, the matter effect is not dominant in the region (A) of 
z/^ channel. 

VII. CP VIOLATION VS. MATTER EFFECT IN LONG-BASELINE NEUTRINO 
OSCILLATION EXPERIMENTS: A NUMERICAL ANALYSIS 

In this section we present the results of our numerical analysis using the two sets of 
parameters (a) and (b) given in (|45|) and (^5]). We do this first for z/^ v^. channel, and 
second for Uj. channel. All the calculations are carried out by using the exact analytic 

expressions found in p8l with the procedure mentioned at the end of section V. 



A. ^ I'e and i/^ channels 

First let us take the parameter set (a). In Fig. 4 we plot -P(z/^ — > z/g) and -P(z/^ — > z/g) 
and the corresponding AP(z/^ z/g) x 100 as a function of ^vn^ / E for different values of s^2- 
We fix the distance L = 250 and 730 km in Figs. 4(a) and 4(b), respectively. We see that 
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from (iv) and (vi) in Fig. 4 that AP(z/^ v^) can be as large as ~ 10 % for L = 250 km 
and ~ 25 % for L = 730 km due to the the matter effect whereas the intrinsic CP violation 
is at most ~ 1.5 %, in agreement with our estimation in the previous section (see Table 1). 
We can conclude that the matter effect dominates over genuine CP violating effect except 
at the exceptional point — As we can see from the second term in eq. (HTl) that 
the matter potential is multiplied by cos2^i2 and hence the matter effect in AP(z/^ — > z/g) 
is suppressed if is close to 0.5. This suppression can also be seen in Fig. 3 in the first 
reference in |T^. We also notice that this factor cos2^i2 gives rise to the sign difference in 
AP(z/^ — > Ue) for > 0-5 and 5^2 < 0.5 as we can confirm from (iv) and (vi) in Fig. 4. 

In Fig. 5 we plot the contour of AP(z/^ — > z/g) on the sl2~ ^f^^ / E plane to see the global 
features of the coexisting CP and the matter effects. We see that the contours for the pure 
CP and the pure matter cases are very different and the amplitude of the latter is larger. 
It may be difficult to see the genuine CP violating effect just by observing AP(z/^ z/g) 
apart from the exceptional region ~ 0.5. 

Let us turn to the parameter set (b). In Fig. 6 we plot as in Fig. 4 P(z/^ z/g) (and 
for v) and AP(z/^ z/g) but for the parameter set (b). In this case the probabilities of 
and z/^ Ve are small, of the order of 1 %, as we can see in Fig. 6 and they may 
be denoted as the minor channels. In this case the matter effect contamination in AP is 
negligibly small, in agreement with our estimation done in the previous section. We confirm 
that the genuine CP violating effect dominates over the matter effect at the parameter (b). 



B. Vt and Vj- channels 

In Fig. 7 we plot P and AP as in Fig. 4 but for the parameter set (a) in the z/^ 
and Vy, —>■ Vr channels. These channels are also minor since P < 2 % but AP^^ is relatively 
large ~ 1 %. We see that the genuine CP violating effect is larger than the matter effect. 
In Fig. 8 we plot as in Fig. 5 the contour of AP(z/^ v^-) on the 5^2 — ^m? /E plane. We 
confirm from these contours that the genuine CP effect is larger than the matter effect, in 
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agreement with our estimation given in Table 2. 

In Fig. 9 we plot P and AP for the parameter set (b). In this case i^^ — > Vr is the 
dominant channel but AP is small ~ 0.7 %. We see that apart from the small s^2 region 
for the basehne L — 730 km the genuine CP violating effect is larger than the matter effect. 
Hence, we conclude that for and i?^ — > channels are relatively free from the 

matter effect. In Fig. 10 we plot the contour of AP(i/^ — > Vt) for this parameter set. We see 
that also for this case the genuine CP effect is larger than the matter effect in agreement 
with our prediction (see Table 2). 

VIII. CONCLUSION 

We have investigated in detail the CP violation in long-baseline neutrino oscillation 
experiments in the presence of matter effect under the assumption of neutrino mass hierarchy, 
Am^ ~ 10~^ eV^ and AM^ ~ a few eV^, motivated by the atmospheric neutrino anomaly 
and the cosmological dark matter. We developed the matter perturbation theory using the 
hierarchy in energy scales, a <^ ^"^^ derived the approximate analytic expressions of 
V — V difference in oscillation probabilities AP(i/^ — > v^). We have found that in a good 
approximation AP can be expressed as a sum of three terms which represent the genuine CP 
violating effect and two different correction terms due to the matter effect. These analytic 
expressions of the AP are useful in understanding the global features of the competition of 
two effects. 

We have studied the question of how large is the magnitude of CP violations due to 
intrinsic CP- violating phase and to the matter effect in the earth. The assumed mass 
hierarchy mentioned above and the interpretation of the atmospheric neutrino anomaly in 
terms of the neutrino oscillations allow us to restrict ourselves into the parameter regions 
strongly constrained by the reactor and the accelerator experiments, (A) small-Si3 and small- 
S23 and (B) Iarge-Si3 and arbitrary S23- 

We have found the following structure (as summarized in Tables 1 and 2). 
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J^fi channel: In the region (A) the matter effect contamination is much larger than 

the genuine CP violation effect except for the case S12 — 0.5. Whereas in (B) the the matter 
effect contamination is negligible compared with genuine CP violation. Unfortunately, the 
magnitude of the latter is small, ~ 0.5 %. 

z/^ — > channel: The matter contamination is smaller than the genuine CP violation 
effect and is not harmful in both regions (A) and (B). The magnitude of the CP violation 
is again not sizable, ~ 1 %. 

Thus, if the mass hierarchy motivated by the dark matter is the truth in nature, it 
appears to be necessary to invent new method for measuring CP violation of ~ 1 % level in 
long-baseline neutrino oscillation experiments. We are planning to discuss an idea toward 



the goal elsewhere |29 



Note added: While we were to complete this paper we became aware of the paper by Arafune 



et al. 10 which addresses the similar topics. However, they consider different neutrino mass 
spectrum and employ a different approximation scheme from ours which requires that < 1 
and -^^^ <^ 1, whereas we only need ^ 1- 

We also note that after this paper had been submitted for publication, a new data on 
atmospheric neutrinos from Super-Kamiokande experiment has appeared |31[]. The data 



seem to favor smaller values of Am^ than the one we assumed in this paper. However, the 
results of the analysis done in this paper is presented so that they are useful for such smaller 
values of Am^. We stress that the perturbative treatment of the matter effect is still valid 
even for Am^ as small as 10~^ eV^. 
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APPENDIX 



We verify that the adiabatic approximation which we have employed to obtain ([T7| ) is in 
fact a very good approximation for the long-baseline neutrino experiments. The adiabaticity 
condition is nothing but the condition 



-^V-^(x) ■ V(x) 
ax 



d)ii\ 



(50) 



To first order in matter perturbation theory the off-diagonal elements of the LHS of (|50[) 
can be expressed as 

" d 



dx 



V^{x) ■ V{x) 



U*-U ■ 

I , s '-^ ej '-^ 



for i ^ j and the diagonal elements vanish for i = j. Then the adiabativicity condition can 
be written as 



\U:,Ue^\ < 



1 



a'{x) 



\ 2E 



(51) 



where the RHS is replaced by Am , the smallest Am^j. The expression ([5T| ) in fact represents 
the sufficient condition, but we will see that it is well satisfied. 

The density gradient in the continental structure is at most 0.2g/cm^ over the depth of 
20 km p6|, which amount to the baseline of 1000 km. Using this a'{x) / a{x) can be estimated 

as 



a'{x) 



a{x) 



1.5 X lO^^km"^ 



-1 



Then, the the adiabaticity condition reads 

IC/- C/.l « 0.83 X 10' (^)' {^f ( ,,,„_.;',„_3,,„- j 

It is clear that the adiabaticity condition is well satisfied in the long baseline experiments 
due to the slow variation of matter density in the continental structure. 
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Region 




r<cp 


/^rnatterl 


/^matter2 


(A) 


~ 1 


~ 10-2 


~ 0.1 


~ 10-2 (0.1) 


(B) 


- 10-2 


- 10-2 


~ 10-4 


~ 10-6 (10-5) 



Table 1: The order of magnitude estimate of maximal values of P, C^^ and (7"*"**^^!' 

Vf, channel in the parameter region (A) and (B) for £" ~ 1 GeV and iSw? ~ 10-2 
For C''^«**er2 g^igQ gjiQ^ ^jig parentheses the values for the case i^vr? ~ 10"^ eV2. 



Region 






r^matterl 




(A) 


~ 10-2 


~ 10-2 


~ 10-3 


~ 10-4 (10-3) 


(B) 


- 1 


- 10-2 


- 10-3 


- 10-4 (10-3) 



Table 2: Same as Table 1 but for channel. 
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4JxlOO 




Fig. 1: Contour plot of 4Jxl00 is shown where J = 012^120235230^3513 sin (5. We set 6^12 = 7r/4 
and 6 = tt/2 which maximize J. The area indicated by dashed hnes with arrows is the region 
excluded by accelerator and reactor experiments at 90 % C. L. obtained in ref. for the 
case AM^ = 5 eV^. The points indicated by the (a) open circle and (b) square are the 
parameters we will use in this work. 



30 



(A-1) 
V3~Vx 



AM^>0 



* » 2 

Am 



(A-2) 



AM^<0 



V3~V, 



(B-1) 

V3~Ve 



AM^> 



^ A 2 



(B-2) 



AM^<0 



Vs-Ve 



Fig. 2: The neutrino mass spectra which we consider in this paper. We assume that 
AM^ = ml — ml c::^ ml — ml ^ a, few eV^ and Am^ = m\ — m\r^ 10~^ eV^. The figure (A- 
1,2) and (B-1, 2) are for the case where mixing parameters are in the region (A) si3,S23 <^ 1 
and in the region (B) S13 ~ 1 and arbitrary S23, respectively. 
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Fig. 3: We plot both the exact and approximated values of AP(t'^ — > i/g) (upper two 
panels) and AP(i/^ — > v^.) (lower two panels) as a function of distance L from the neutrino 
source. We fixed the mixing parameters as s\.^ = 3.0 x 10~^, = 2.0 x 10~^ for the left 
two panels (a) and (c) and = 2.0 x 10~^, 5^3 = 0.98 for the right two panels (b) and 
(d). The remaining parameters are fixed to be the same for all the case (a-d), i.e., — 0-3, 
AM^ — 5 eV^ and 5 = 7r/2. The solid lines and open circles are for the case l^rri^ /E — 10~^ 
eV^/GeV and the dotted and open diamonds are for the case iSrv? jE = 10~^ eV^/ GeV. 
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Fig. 4: We plot in (i), (ii) and (iii) -P(t'yn — > fg) and -P(P^ — > Pg) for the parameter set 
(a) = 3.0 X 10-3, s?3 = 2.0 x 10-^ as a function of Am'^/E for sj^ = 0.2, 0.5 and 0.8, 
respectively. We fixed the remaining parameters as AM^ = 5 eV^, 5 — 7r/2. The figures (a) 
and (b) are for L — 250 km (a) and L — 730 km (b), respectively. In (iv), (v) and (vi) we 
plot the corresponding AP(i/^ v^) x 100 and the different curves are for the cases where 
only the CP (dotted line), only the matter (dashed line), and both the CP and the matter 
(solid fine) effects are taken into account. 
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Fig. 5: We plot, for the parameter set (a) s^g = 3.0 x 10"^ and sfg = 2.0 x 10~^, the contour 
of AP(z/^ — > z/e) X 100 in the — ^rn? / E plane for the cases where only the CP effect (a, 
d), only the matter effect (b, e), and both the CP and the matter effects (c, f) are taken into 
account. The upper and the lower three panels are for L = 250 km and 730 km, respectively. 
In (c, f) we also indicate by the long-dashed symbol the line = 0.5 where the matter 
effect vanishes. The remaining parameters are fixed to be the same as in Fig. 4. 
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(a) L = 250 km 
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Fig. 6: Same as in Fig. 4 but for the parameter set (b) S23 — 2-0 x 10 ^, 5^3 = 0.98. 
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(b) L = 730 km 
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Fig. 7: Same as in Fig. 4 with the parameter set (a) S23 = 3.0 x 10 ^, 5^3 = 2.0 x 10 
for i/„ — > Ut and Pu — > channels. 
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Fig. 8: Same as Fig. 5 but for i/^ — > Ut channel with the parameter set (a) S23 — 3-0 x 10 , 



sfg = 2.0 X 10-^ 
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(a) L = 250 km 
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Fig. 9: Same as in Fig. 7 but for the parameter set (b) S23 — 2-0 x 10 ^, 5^3 = 0.98. 
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(d) Only CP (L=730 km) 
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Fig. 10: Same as Fig. 8 but for the parameter set (b) sig = 2.0 x 10 ^, sj^ = 0.98. 
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